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We present the calculations of the Curie temperature and magnetization of doped EuO both in 
the absence and in the presence of external magnetic field. The calculations were performed both 
for the free electrons model and for the model with finite electron band width. Both models give 
similar results for the magnetization, close to Brillouin function. 



I. INTRODUCTION 

The perspectives of the development of spintronics has 
led to renewed interest in rear-earth oxide ferromagnetic 
semiconductors such as europium chalcogenides, in par- 
ticular EuO and EuS, which have a rock salt (fee) struc- 
ture, with a lattice constant of 5.14 A for EuO and 5.96 
A for EuS, and whose magnetism arises from partially 
filled and highly localized 4/ states. Stoichiometric EuO 
and EuS are regarded as typical Heisenberg ferromag- 
nets, with a Curie temperature (X"c) 69.8 K for EuO 
and 16.6 K for EuS respectively [TJ |2J. From the point of 
view of electronic properties these materials are semicon- 
ductors, with a band gap at room temperature of 1.12 
eV for EuO and .6 eV for EuS. For applications EuO 
looks especially promising [3J; it has the third strongest 
saturation magnetization of all known ferromagnets [I], 
one of the largest magneto-optic Kerr effects [5], pro- 
nounced insulator-to-metal transitions [SHI] as well as 
colossal magnetoresistance effects [10]. For the use of 
these materials in spintronics applications it is desirable 
to increase Tq as much as possible. This can be achieved, 
in particular, by doping EuO with the rear earth metals 
with valency 3 (Gd [3J, La |llj ) or with oxygen vacan- 
cies. Introduction of carriers into the conduction band 
by doping leads to indirect exchange between the local- 
ized spins, thus the Curie temperature depends upon the 
electron concentration. 

The seminal calculations of the influence of the indirect 
exchange on magnetic properties of doped Europium ox- 
ides were performed by Mauger [12]. In these calcula- 
tions the s — f interaction is treated in the second order 
of perturbation theory, very much similar to the RKKY 
theory [21 . However, the RKKY approach was substan- 
tially modified due to some specific features of magnetic 
semiconductors. First, in magnetic semiconductors the 
electronic gas is not always degenerate, contrary to met- 
als. Second, due to the fact that the atomic exchange 
is not small compared to Fermi energy, the temperature 
dependence of the indirect exchange cannot be neglected. 
Thus the calculations were performed explicitly at a finite 
temperature T and not at T = as usual. Also the con- 
duction band of a finite width was considered, contrary 
to the RKKY assumption of the free electrons dispersion 
law. Mauger calculations still remain an important refer- 



ence point for the experimentalists in the field, including 
recent experiments [52] . 

The theory of indirect exchange was substantially 
modified by Nolting [13] [H] [15] . The modified the- 
ory involves 3 types of correlation functions, describ- 
ing itinerant electrons correlations, local-spin correla- 
tions and mixed itinerant-electron-local-momcnt corre- 
lations. This theory is able to describe magnetic semi- 
conductors (EuO, EuS), diluted magnetic semiconduc- 
tors (Gai_ x Mn x As), magnetic metals (Gd, Dy,Tb) and 
CMR ma,teiia\s(Lai- x Ca x Mn03). Using this model 
Nolting had performed calculations of density of states 
(DOS) and energy band structure for EuO in T = and 
finite temperature. |16j |17) 

However, in the present work we decided to limit our- 
selves with the Mauger type calculations. Our aim was 
to check up to what extent the final results arc influenced 
by the dispersion law of the itinerant electrons. We com- 
pared the results obtained in the framework of the free 
electrons model with those obtained in the framework 
of the model with the finite electron band width used 
by Mauger [T5]. As an additional modification we per- 
formed the calculations using Matsubara Green functions 
to calculate indirect exchange at a finite temperature. 



II. INDIRECT EXCHANGE BY MATSUBARA 
GREEN FUNCTIONS 

The Hamiltonian we start from is 
H = -^ n Y l f ^(r)VVa(r)d 3 r-^/(R ij )S i S j 

a ij 

+ J d/ ^S I a(R l ), (1) 

i 

where a(r) = iil(r)<T al3 ip a (r). 

We'll use relevant for EuO approximation of weak s — f 
exchange coupling, and calculate the indirect exchange 
between localized spins appearing due to their inter- 
action with conduction electrons in the leading order 
of perturbation theory. (For the intermediate coupling 
regime see e.g. O US].) Though our aim is to obtain 
finite temperature results for the indirect exchange, we'll 
start from the derivation of the RKKY interaction for 
T = as presented in Ref. [55]. We can write down spin 
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polarization in the Fermi gas c l (r) = (ipa{ T ) a 'ap' l Pp( T ) 
using Green function 



G(iu n ,r) = 



d 3 p 



(2?r) 3 iuj n - C P + ^sign w, 



, (11) 



lim (-ia- % a gGp a (T, t; r', t')) , 

f'->t+0,r=r' v p ' 

(2) 



where the summation on indices a, ft is implied. Let 
assume that there is a localized (at r = 0) spin S 
interacting with the local spin density of conduction 
electrons 



H int = J df S l a l (r = 0), 



(3) 



where a 1 = ^ a (r)a % a ^i\) ^{v) . l n tne nrs l or der with 
respect to interaction 

G (1) (e,r),r') = J d/ 5V^G (e,r)G (e, -r'). (4) 



Hence 



<x*(r) = -2U df S l / G" (e,r) 



de 
'2tt' 



(5) 



where £ p = #- — -Ep and oj„ = (2n + 1)ttT. (We have 
taken into account scattering of conduction electrons 
[27] : t is the scattering time.) Making fraction decom- 
position we obtain 



G(iu n ,r) 



47T 2 7 



1 



1 



K — p K + p 



sin pr dp 



2irr 

where 

k = (/2m(i?f + ioj. 



(12) 



^rsign u) n ). (We use the value of 
the root which has the same sign of the imaginary part 
as U3 n .) Finally we obtain 



a\r) = 2J rf/ ^T^G 2 ( 



(13) 



Eq. |13| leads to 

Jeff (Rij) = 2J 2 df Tj2G 2 (^n, Rii). (14) 



III. FREE ELECTRONS MODEL 



For the free electrons model we can take 



G (e,p) = — ^ — ITT^ ' ( 6 ) 

e + tip — p z / Zm + iosigne 



and hence 



Go(e,r) 



m 
2nr 



^isigne^r 



where k — ^/2m(E F + e + i&igne) and calculate the 
integral ([5| to obtain 

,., . :2mk% /cos x sinx\ . . 

a\r) = J d fS*—^ ( — ), (8) 



where x = 2k F r. Thus we obtain the well known RKKY 
result 



where 



H 'ex - - y ^2 J eff^ij) S i S h 



23% {k F a /2f 

J e ff(Kij — 5 

fj, 7T 3 

ain(2k F Rij) — 2k F Rij cos(2k F Rij) 
(2k F R l3 ) i • 



(9) 



(10) 



where /i = k F /2m and A = kpTjm* . To make calcula- 
tions for finite temperature [25] , first we have to calculate 



To take into account spin polarization of the conduction 
electrons we should modify Eq. (fl4| to 



Jeff(Rij) - Jdf T J! G a( iuJ n,R 



IJJ) 



where 



and k± 



G a {iw n ,r) 



2irr 



(15) 



(16) 



2m(E F ± JdfSa/2 



^ sign u>„ 



(7) 

where a =< S z > /S is the reduced magnetization. 



IV. FINITE ELECTRON BAND WIDTH MODEL 



In this model the dispersion law is given by: 



W 

Ek = -y (1 _ cos ( ka )) 



(17) 



where a is the lattice constant ,W is the conduction 
bandwidth and k is the wave vector. The effective 



exchange in Eq.(15) is replaced by: 

2 



Jeff(Rij) 



+.f (E k - J d fSa/2) 



ksin (kRij) k' sin (k'Rij) 
cos (k'a) — cos (ka) 



dkdk' 



(18) 



where P means the Cauchy principal part of the integral 
over k, z is the degeneracy of the conduction band (as- 
sumed to be 2 in our case), not including the spin. V is 



3 



the volume of the crystal and / (E) is Fermi function. 



EFFECTIVE SPIN HAMILTONIAN 



Following Mauger we chose Jdf — 0.13 eV for EuO. We 
take density of states effective mass m being equal to the 
free electron mass 129 



The calculations, the results of which are presented 
below, are based on the effective Hamiltonian 

H = /(R^SiS; - i£ Je//(Ry)S?S|,(19) 



where J e ff(Rij) is given by Eq.(15l or by Eq.(18|. 
As one can see, for the case of polarized electrons 
we are getting some anisotropic Heisenberg model for 
the indirect exchange. However, in the mean field 
approximation(MFA) only z direction is relevant. The 



sum in Eq.(15l can be calculated only numerically. 
Further on we'll present the calculations corresponding 
to r = oo, leaving the analysis of the influence of the 
finite scattering times for the future publication [2"B] . 
The emergent magnetic problem is solved in the mean 
field approximation. The equation for the reduced 
magnetization a is [T2"] 



a = B, 



aS 2 (1(0) + £\ J e //(%) J + gfi B SB, 
T 



(20) 



where Bg is the Brillouin function, g is the g factor, [Lb is 
Bohr magneton and B ext is the external magnetic field. 
The quantity 1(0) we chose to fit the Tq of the stoichio- 
metric compound. Because the average magnetization of 
localized spins enters into the integral, the calculations 
should be done selfconsistently. We are taking in consid- 
eration the Eu ions have a structure form of FCC (Face 
Center Cubic) lattice and the sum over i£y is done ac- 
cording to Tabic I. 
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24 
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TABLE I. FCC Number of neighbors and Distance between 
neighbors. 



VI. RESULTS 



A. Free electrons model 



First, we calculated the magnetization a as a function of 
temperature for different electron concentrations accord- 



ing to Eq.(20) in the absence of external magnetic field. 
The results are presented in Fig(T]and[2] 
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FIG. 1. The magnetization as a function of temperature (T) 
for Jdf = 0.13 eV and for different electron concentrations: 
n = (Thick line), n = 0.05 (Dashed line), n = 0.1 (Dotted 
line) and n = 0.15 (DotDashed line). 




FIG. 2. The magnetization as a function of temperature (T) 
for Jdf = 0.13 eV and for different electron concentrations: 
n = 0.2 (Thick line), n = 0.25 (Dashed line), n = 0.3 (Dotted 
line) and n = 0.35 (DotDashed line). 



FigjT] and [2] shows clearly that the magnetization is like 
Brillouin function even for n ^ 0. we see also how the in- 
direct exchange (n ^ 0) changes the critical temperature 
compared to Curie temperature (69.8K) of stoichiomet- 
ric compound (n = 0). We plotted Curie temperature 
Tq as a function of electron concentration n on Fig(3j 
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FIG. 3. Curie temperature (lb) as a function of electron 
concentration (n) for Jdf = 0.13 eV where the points represent 
real points and the line is interpolation. 



FIG. 5. The magnetization as a function of temperature (T) 
for Jdf = 0.13 eV and n — 0.1 and for different external mag- 
netic fields: B ext = (Thick line), B ex ± = 1 Tesla (Dashed 
line), B ext = 5 Tesla (Dotted line) and B ext = 10 Tesla (Dot- 
Dashed line). 



In Fig|3] we see that Curie temperature has maxi- 
mum value (1A9K) at electron concentration of n = 
0.05 and return to the temperature of stoichiometric 
compound(69.8i ; ir ) at electron concentration of n — 0.3. 
For n = 0.45 the Curie temperature drops to zero, which 
means that above this concentration the ferromagnetic 
phase does not exist even at T = 0. 
The reason for the suppression in Curie temperature for 
large electron concentrations is that as the doping or oxy- 
gen vacancies level increases, the wavelength of the effec- 
tive exchange oscillations (1/kp) becomes shorter and, 
hence, increasingly, anti-ferromagnetic, which ultimately 
suppresses the ferromagnetic transition. 



Now we use Eq. ( 20 ) to calculate the magnetization a as 



a function of temperature for different electron concen- 
trations and for different external magnetic fields. The 
results are presented in Fig. Hp] 
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FIG. 6. The magnetization as a function of temperature (T) 
for Jdf = 0.13 eV and n = 0.3 and for different external mag- 
netic fields: B ext = (Thick line), B ex t = 1 Tesla (Dashed 
line), B ext = 5 Tesla (Dotted line) and B ex t = 10 Tesla (Dot- 
Dashed line). 
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FIG. 4. The magnetization as a function of temperature (T) 
for Jdf = 0.13 eV and n — and for different external mag- 
netic fields: B ext = (Thick line), B ex t = 1 Tesla (Dashed 
line), B ext = 5 Tesla (Dotted line) and B ex t = 10 Tesla (Dot- 
Dashed line). 



Fig|4][6] shows clearly that the magnetization is not like 
brillouin function for « / because of the existence of 
an external magnetic field. 



B. Finite electron band width model 

First, as before, we calculated the magnetization a as 
a function of temperature for different electron concen- 



trations according to Eq.(20) in the absence of external 



magnetic field. The results are presented in Fig(7]and[8] 
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FIG. 7. The magnetization as a function of temperature (T) 
according to free electrons model for Jdf = 0.13 eV and for 
different electron concentrations: n — (Black Thick line), 
n = 0.005 (Black Dashed line), n = 0.015 (Black Dotted 
line), n = 0.028 (Black DotDashed line), n = 0.044 (Gray 
Thick line) and n — 0.065 (Gray Dotted line). 
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imum value (155 K) around electron concentration of 
n = 0.1 and return to the temperature of stoichiomet- 
ric compound (69.8K) around electron concentration of 
n = 0.25. 

Like before, the reason for the suppression in Curie 
temperature for large electron concentrations is that 
as the doping or oxygen vacancies level increases, 
the wavelength of the effective exchange oscillations 
(1/kp) becomes shorter and, hence, increasingly, anti- 
ferromagnetic, which ultimately suppresses the ferromag- 
netic transition. 

If we do a comparison between the two models we can see 
that the concentration dependence of the Curie tempera- 
ture is similar: the maximum Curie temperature achieved 
for the optimal doping is almost the same, the concen- 
trations of electrons corresponding to optimal doping are 
very close. However, for the free electrons model (Fig. [3]) 
the decrease of the Curie temperature with the increase 
of electron concentration is slower than in the finite elec- 
tron band width model (Fig. [9]). 

As before, we calculated the magnetization a as a func- 
tion of temperature for different electron concentrations 
and for different external magnetic fields. The results are 
presented in Fig. [lO]fT2 
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FIG. 8. The magnetization as a function of temperature (T) 
according to Mauger model for Jdf = 0.13 eV and for different 
electron concentrations: n = 0.091 (Black Thick line), n = 
0.11 (Black Dashed line), n = 0.14 (Black Dotted line), n = 
0.17 (Black DotDashed line), n = 0.2 (Gray Thick line), n = 
0.25 (Gray Dotted line) and n = 0.29 (Gray Dashed line). 



We plotted Curie temperature Tq as a function of elec- 
tron concentration n on Fig|9] 
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FIG. 9. Curie temperature (Tc) as a function of electron 
concentration (n)- Mauger's result 

In Fig. [9] we see that Curie temperature has max- 
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FIG. 10. The magnetization as a function of temperature 
(T) for J d f = 0.13 eV and n = 0.044 and for different exter- 
nal magnetic fields: B ex t = (Thick line), B ex t = 1 Tesla 
(Dashed line), B ext = 5 Tesla (Dotted line) and B ext = 10 
Tesla (DotDashed line). 



G 




60 120 180 240 300 360 



T(°K) 

FIG. 11. The magnetization as a function of temperature 
(T) for Jdf = 0.13 eV and n = 0.14 and for different exter- 
nal magnetic fields: B ex t = (Thick line), B ex t = 1 Tesla 
(Dashed line), B ext = 5 Tesla (Dotted line) and B ex t = 10 
Tesla (DotDashed line). 
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FIG. 12. The magnetization as a function of temperature 
(T) for Jdf = 0.13 eV and n = 0.2 and for different exter- 
nal magnetic fields: B ex t = (Thick line), B ex t = 1 Tesla 
(Dashed line), B ext = 5 Tesla (Dotted line) and B ex t = 10 
Tesla (DotDashed line). 

Again, we can see clearly that the magnetization is not 
like Brillouin function for n ^ because of the existence 
of an external magnetic field. 

VII. CONCLUSIONS 

We presented the calculations for the influence of the in- 
direct exchange on the magnetic properties of EuO in 
the framework of Matsubara Green functions technique. 
We used different electron dispersion laws: free electrons 
and finite conduction electron band width. Qualitatively 
the results for the concentration dependence of the Curie 
temperature and temperature dependence of the magne- 
tization were similar, the maximum Curie temperature 
achieved for the optimal doping is the same, the con- 
centrations of electrons corresponding to optimal doping 
are very close. However for the free electrons model the 
decrease of the Curie temperature with the increase of 
electron concentration is slower than in the finite elec- 
tron band width model. The temperature dependence 



of the magnetization in the former case looks closer to 
the Brillouin function. However, all experimental results 
on doped EuO shows a non-like Brillouin magnetization 
curve and it seems to have a second dome in the magne- 
tization curve. Some authors argue [5] [IB] [T^j .therefor, 
that the second dome represents a second critical tem- 
perature {Tc). However, in our calculations we did not 
see any deviation from Brillouin function for doped EuO 
in the absence of an external magnetic field .One of the 
possible explanations for this contradiction could be the 
existence of an external magnetic field in experimental 
systems. 

Also, some experiments show values of Curie tempera- 
ture which are larger than our results {Tc = 170K [18] 
and Tc = 180K )■ Nevertheless , we emphasize that 
both models we used , despite their simplicity, have a good 
match to most of experimental results of temperature vs. 
electron concentration curve and the high Curie temper- 
ature values we have mentioned, are exceptional . 
In this paper, we explained the existence of maximum 
of Tc by ferromagnetic and anti-ferromagnetic ordering 
which increases and reduces respectively the indirect ex- 
change. However, some experimental works [3] show the 
fact that only a small fraction of the introduced dopants 
act as a donor. These results could provide an alterna- 
tive explanation the for saturation (or maximum) of Tc 
which stems from the saturation and decreasing of the 
concentration of "mobile" charge carriers. 
The problem of the influence of doping was recently an- 
alyzed in Ref. 30 . It was shown there that there are 
two competing factors. On one hand, like it was tradi- 
tionally considered, free carriers induce the Ruderman- 
Kittel-Kasuya-Yosida (RKKY) interaction. On the other 
hand, since the bottom of the conduction band consists 
mainly of the majority spin, the doped electron will en- 
ter the spin-polarized manifold, and this results in the 
onset of moment in the 5d band. Both these factors were 
studied in Ref. [3U] for EuO using the virtual crystal 
approximation. The second factor was not taken into ac- 
count in the main body of our paper (and the first was 
taken into account in a mode advanced approximation, 
than virtual crystal). 

The comparison of the results of two approaches (free 
electron model and finite electron band width model) 
between themselves and with the experiment can give 
an idea, what features in the observed experimental be- 
havior are robust (the dispersion law independent), and 
what aren't. Thus it may help to understand the physics 
of magnetic semiconductors. 
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